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1. Introduction 

In this work we want to test the matching between Lattice QCD and ChPT in a special finite 
volume regime in which the valence quarks have masses tn v in the £ -regime (that is wt v £V ^ 1, but 
LA X 1 where A x is the cutoff of ChPT, £ is the quark condensate, V = L 3 x T is the volume). 
As discussed in [1], under these conditions slow pion modes dominate the path integral in the 
effective theory, so the expansion in powers of M^/A 2 ~ p 2 /A 2 breaks down and a new power- 
counting has to be used. In particular the pion zero modes have to be treated non perturbatively. 
Nevertheless, the low energy constants (LECs) that appear in the ChPT Lagrangian are unaffected 
by finite size effects so that this regime provides a framework in which the LECs can be determined 
with different systematics with respect to usual "infinite volume" studies. 

Our aim is to use a mixed action approach, in which chiral symmetry is exactly preserved at the 
level of valence quarks only. We consider overlap valence quarks on top of Nf = 2 CLS ensembles, 1 
obtained from simulations with non-perturbatively 0(ct)-improved Wilson sea quarks. The use 
of cheaper Wilson quarks allows us to consider rather small lattice spacings (a ~ 0.08/m) and 
sufficiently large physical volumes (L ~ 2fm). We also take a Partially Quenched approximation 
in which the sea quark masses m s stay in the /^-regime, that is they satisfy m s LV S> 1 - A consistent 
power counting for this regime was introduced in [2]. 

The use of overlap in the valence sector is clearly suitable to study low energy QCD phe- 
nomenology, where the effects of spontaneous breaking of chiral symmetry are more visible. More- 
over it allows us to define unambiguosly the topological charge v through the index theorem. As a 
first application we have matched the dependence of (v 2 ) on the sea quark mass to the prediction 
of NLO ChPT in the mixed regime, that depends on the LECS: £ and a combination of Lg, L-j and 

In the £ -regime at a fixed v, the LO partition function of ChPT has been shown to be equiva- 
lent to a random matrix theory (RMT) [3], where many analytical predictions can be obtained for 
spectral quantities, such as the spectral density or the distribution of individual eigenvalues of the 
Dirac operator. Here we show that ChPT in the mixed regime, up to NLO, can also be matched to 
a RMT. The dependence on the sea quark mass of the parameter of the RMT can be extracted from 
this matching and can be shown to depend on £ and L§. 

A detailed paper has been published recently [4], to which we refer for a more extensive 
discussion. 



2. Matching QCD in the mixed-regime with Random Matrix Theory 

As argued in [2], the heavier ^-regime sea quarks behave as decoupling particles that can be 
integrated out. The dependence on m s can then be reabsorbed in the LECs. The matching to a RMT 
requires the further integration out of the non-zero momentum modes in ChPT. 

For a theory with Nf quarks in the £ -regime the partition function of ChPT at fixed V, iF v c/iPr 
can be written at LO as: 



yChPT 
'V 



MQ= \d%\ [^o]det(^o) v expf^Tr[^(^o + ^o t )lV _J ^ YTr[ ' 9 ''^^ ] W 
J JU(Nf) \ 2 L J J 



'https://twiki.cern.ch/twiki/bin/view/CLSAVebHome 
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where Ot = 5,-yra, is the Nf x Nf mass matrix, t/o and E, parametrize respectively the non perturba- 
tive zero modes and the perturbative non zero modes of the pions. In finite volume the momentum 
is quantized and for sufficiently low energies the non zero modes can be decoupled, order by order 
in the £-expansion if the Pseudo-Goldstone masses M 2 ab = £(m ^ m " ) satisfy M ab «r'. This leads 
to the definition of a Zero Modes Chiral Theory (ZMChT) with partition functional 

3^ chT (JK) « / [dUo] det(U ) v exp ( ^Tr \jt(U + t/ ( j)l ) . (2.2) 

In [3], it was shown that !^™ chl (at) is equivalent to the partition functional of a Random Matrix 
Theory (RMT) of matrices W of size N, that depends only on the number of flavours, Nf, and the 
corresponding mass parameters m, viz. 

^™ chT (at) ~ RMT N {N f ,m} , (2.3) 

and there must be an identification Nriii = miLV = From this relation spectral correlation func- 
tions of the Dirac operator can be related to the corresponding quantities computed in RMT. In 
particular, if x is the k-th smallest eigenvalue of the matrix v W^W, the probability distribution 
associated to the microscopic eigenvalue £ = Nx, pY(^;{^l}) can be calculated in RMT. Then also 
the distribution of individual low-lying eigenvalues of the massless Dirac operator can be predicted 
[5] from the equivalence Eq. (2.3) in terms of just one parameter, the chiral condensate E. 

This relation has been tested in the quenched approximation [6], respectively in Nf = 2 and 
Nf = 2 + 1 dynamical simulations in [7] and [8]. 

In the mixed regime with N v £ -regime quarks (M vv <C L _1 ) and N s /^-regime quarks (L _1 ~ 
M ss ), at low energies all non-zero modes and also the zero modes of the heavier pions must be 
integrated out. A ZMChT is obtained corresponding to Nf = N v flavours as opposed to the total 
N v +N s that one would obtain if all the quarks were in the £-regime. 

At LO the ZMChT in the mixed-regime is simply ChPT without the heavy modes (the inte- 
gration over them gives an irrelevant normalization factor): 

Z ™m hT \w x j UN ) [^o](det[/ ) v exp (^Tr [at y (t/o + f/j)]) ■ ( 2 - 4 ) 

where is the mass matrix at projected to the valence sector. According to the eq. (2.3), this 
partition function is then equivalent to an ./V v RMT. In particular, it is important to stress that the 
ZMChT has 7Y V flavours, while the full ChPT from which it is derived corresponds to Nf = N v +N S 
flavours. In particular, for Af v — > 0, the ZMChT or RMT we expect to find is the quenched one, 
while the couplings should be those of an Nf = N„ theory. 

The matching at NLO still does not modify the structure of the ZMChT theory. The dynamics 
of heavy modes can be absorbed in the couplings appearing at LO in eq. (2.4), that is £. The 
integrations over the perturbative modes result in a change £ — > E e ff, which for degenerate quarks 
and in the limit A^ v — > 0, reads: 



I 2mJL 
E eff (m s ) =Xi 1 + -^- 



T + wW 1/4 ) +3 ^W-^(]^) 



Pi 

2F 2 V 1 / 2 




(2.5) 
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where M p = 770 MeV is the p meson mass, while the function gi (M, L, T) expressing the volume 
effects on the propagator and the /3, are defined in [9]. In our case T/L = 2, j3i = 0.08360 and 
[$2 = —0.01295. The quenched limit of the zero-mode integral over U (N v ) is calculated in [10] and 
matches quenched RMT (qRMT). 

Comparing the eigenvalues computed numerically in this Partially Quenched setup with the pre- 
dictions of qRMT, we can extract E e g of eq. (2.5): 

(^qRMT = ^WImVW^W, (Ck)qRMT = / CpW'M- ^ 

On the other hand, the prediction for the ratios (C/Oormt/ (^)qRMT * s parameter-free and can be 
compared directly with (A^)q CD / (A/)q CD at any fixed M ss . 

The distribution of topological charge is, on average, controlled by the sea quarks only. In the 
case when N v — > and the N s sea quarks are degenerate, we found: 



2 m s iy 
(v ) = 



Nf — 1 / Lm s ( 2m JL \ / 2m JL 



32£»7 

+ —^(Ll(M p )+N s L r 6 (M p )+N s L r 7 (M p )) 



(2.7) 



This result agrees with [11] in infinite volume and with [12] in the ^-regime. Note the appearance 
of L' 7 (M P ), for which no prediction has yet been obtained on the lattice. 



3. Results on Dirac spectral observables 

We have carried out our computations on the CLS lattices of size 48 x 24 3 named D4, D5 
and Dg[13]. The configurations have been generated with non-perturbatively 0(a) -improved Wil- 
son fermions at /3 =5.3 using the DD-HMC algorithm and saved configurations are separated 
by 30 HMC trajectories of length % = 0.5. The lattice spacing has been determined to be a = 
0.0784(10) fm in [13], which implies that our lattices have physical size L ~ 1.88 fm and sea pion 
masses of 426, 377 and 297 MeV, respectively. However, preliminary results from other determi- 
nations (eg using the mass of the CO baryon) yield a ~ 0.070 fm [14]. We will consider both values 
in our analysis. 

On these configurations we have built the massless Neuberger-Dirac operator Dn [15] with 
the parameter s that governs the locality fixed to s = 0.4. We have computed v using the index 
theorem v = n + — n_ where n + («_) is the number of zero modes of Dn with positive (negative) 
chirality, and the 10 lowest eigenvalues 7 of Dn- Given that the 7 appear in general in complex 
conjugated pairs and lie on a circle in the complex plane 7 = = ( 1 — e"^ ) , we have used the projection 
X = yjyy* = = y/2(\ — cos0) to compare with RMT. Expectation values were computed at fixed 
absolute value of the topological charge |v|. 

In Fig. 1 (left) we report the ratios (A*) v / (A/) v normalized to the corresponding qRMT pre- 
dictions, for |v| = 0, 1,2 and for several combinations k,l given at the bottom of the plot. This 
allows to appreciate the precision and level of agreement with qRMT of each specific case. While 
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Figure 1: Left: eigenvalue ratios at fixed |v|, normalized to qRMT predictions, for the indices (k,l) 
(2, 1), (3, 1), (4, 1), (3,2), (4,2), (4,3). Right: £ e ff bare extracted from the matching in eq. (3.1), for k 
1,2,3,4 and | v | =0,1,2. The data for D4 have been shifted in the horizontal axis for better clarity. 



the RMT prediction seems to work well for ratios not involving X\ , the ratios (A*) / (Ai ) exhibit 
somewhat more significant deviations. 

In the spirit of the mixed-regime ChPT analysis, our data also allow to study the mass de- 
pendence of the effective condensate, cf. eq. (2.5). The values of the bare effective condensate 
extracted from the matching 

E eff (M,) = i|^ (3.1) 

for k = 1,2,3,4 and |v| = 0, 1,2 are shown in Fig. 1 (right), where one can observe that, at fixed 
value of the sea quark mass, E e g does not depend on k and v within the statistical precision (with 
larger errors for k = 1). We average over k = 2,3,4 and |v| = 0, 1,2 and estimate systematic effects 
by adding the values for k = 1 in the average. 

4. Fits to NLO Chiral Perturbation Theory 

On the basis of the evidence presented in the previous section, now we assume that the match- 
ing to ChPT in the mixed-regime works in the range of masses and volumes we used, and try to 
extract the low-energy couplings from the sea-quark mass dependence of the two quantities £ e g 
and (v 2 ). 

We first consider the topological charge distribution. The statistical error in this quantity is 
fairly large, but it is encouraging to see that there is a very clear dependence on the sea quark mass 
(see Fig. 2 left). We have fitted both to the full NLO formula in eq. (2.7), and to the linear LO 
behaviour. In either case m s is taken to be the PCAC Wilson mass renormalized in the MS scheme 
at 2 GeV. The results for D4 and D5 are taken from [16], while we have computed that of D6- 

At NLO we fit for £ and the combination [Lg + 2 (Lg + Ly ) ] (M p ) . The value of F is taken to be 
90 MeV; the systematic uncertainty related to this choice is estimated by varying F by ±10 MeV. 
In physical units we get for a = 0.078 — 0.070 fm: 

3 



I MS (2 GeV) 



262 < - 33 ^ MeV 
zoz (34)(5) iViev 



287( 35 )( 5 ) MeV 

ZO/ (36)(7) iVleV 



(4.1) 



[Ll + 2(L r 6 +Uj))(M p ) = 0.0018(30) - 0.0023(43) (4.2) 
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Figure 2: Left: (v 2 ) versus the sea quark mass am^ s (2GeV). The smaller errors are the statistical ones and 
the largest include our estimate of autocorrelations. The dashed and solid lines correspond to the best LO 
and NLO ChPT fits respectively (for a = 0.078 fm). Right: £ e ff<2 3 - The dashed and solid lines correspond to 
the best fit for £ and L\{M p ) at F = 90 MeV, taking the scale to be 0.070 fm and 0.078 fm respectively. 



where the first error is coming out of the fit and the second is the effect of changing F. Since our 
ensembles are not large enough to allow a precise estimate of the effect of autocorrelations, we 
quote errors computed with a conservative approach [4]. 

Let us now turn to £ e ff. In this case, the dependence on m s is expected starting at NLO in 
ChPT. We perform a two-parameter NLO fit, where we fix F and fit for E and L r 6 (M p ). 
In order to obtain the scalar density renormalization factor in the MS scheme for the valence overlap 
fermions we apply the matching condition 



7 MS 



-1 



111 



overlap 



Aft 



m MS {2 GeV) 



Wilson 



Aft 1 



(4.3) 



at aM ref = aM ss . We obtain Zf s (2 GeV) = 1.84(10), where the error is dominated by the one in 
aM vv , i.e. in the determination of the unitary point. Obviously, several checks need to be done to 
ensure that this result is robust, such as checking the dependence on the reference pion mass, as 
well as on the sea quark mass. 

With this estimate for the renormalisation factor and fixing F = 90 MeV, the result we obtain 
from the fit is, for a = 0.078 - 0.070 fm 



L M *(2 GeV) 
L r 6 (M p ) 



9 «(12)(1) 
^ JJ (13)(4) 

0.0010(6) 



MeV 



280 (13)(4) MeV 



0.0015(10), 



(4.4) 
(4.5) 



where the only systematic error that has been estimated is that associated to the change of F by 
±10 MeV. 

Although the fit is good (see Fig. 2 right), it would be desirable to have more sea quark masses and 
smaller statistical errors to assess the systematics. Particularly useful would be to test the finite-size 
scaling. 

This value of £ is consistent with the one obtained from the topological susceptibility above, 
and both are in nice agreement with the alternative determination of [17], that extracted the con- 
densate from the spectrum of the Wilson-Dirac operator on Nf = 2 CLS configurations at the same 
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lattice spacing and sea quark masses, but in a larger physical volume. 
5. Conclusions 

Using a mixed action approach, together with being in the mixed regime and in the PQ ap- 
proximation has allowed us to match the Dirac low spectrum with the ChPT predictions with sig- 
nificantly finer lattices than in previous studies. Our determination of £ lies in the same ballpark 
as many of these previous determinations. More work is however required to quantify the system- 
atic uncertainties involved in this method: autocorrelations, systematics in the chiral fit and finite a. 
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